Nonlinear charge injection in organic field-effect transistors 
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Transport properties of a series of poly(3-hexylthiophene) organic field effect transistors with Cr, 
Cu and Au source/drain electrodes were examined over a broad temperature range. The current- 
voltage characteristics of the injecting contacts are extracted from the dependence of conductance 
on channel length. With reasonable parameters, a model of hopping injection into a disordered 
density of localized states, with emphasis on the primary injection event, agrees well with the field 
and the temperature dependence of the data over a broad range of temperatures and gate voltages. 
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I. INTRODUCTION 

Field effect transistors based on organic semiconduc- 
tors (OFETs) have attracted interest for their poten- 
tial applications as inexpensive and flexible electronics 
The physics of charge injection at the metal-organic semi- 
conductor (OSC) contact in such field-effect devices re- 
mains poorly understood. Contact resistance is often 
neglected when inferring the mobility from transistor 
transport characteristics. It is known, however, that 
the parasitic series resistance, R s , at the OSC-metal in- 
terface can play an important roleQ. Several experi- 
ments have already shown that can easily dominate the 
intrinsic channel resistance, Rch, in short channel (few 
microns and below) devices 0,0, 0. Approaches used 
to differentiate between contact and channel resistances 
include analyses of single device characteristics [(J [jj, 
scanning potentiometry]^, 0, 0, and scaling of to- 
tal device resistance with channel length in a series of 
devices H S S El [H [3 . 

In a previous study fU, we reported measurements of 
R s and the true channel resistance in bottom contact 
poly (3-hexylthiophene) (P3HT) field effect transistors 
as a function of temperature. In that study, which used 
Au source and the drain electrodes, we found that R s 
correlates inversely with mobility over four decades, over 
a broad range of temperatures and gate voltages. This is 
consistent with the predictions of a recent theory [T^LIT^ 
of OSC-metal contacts incorporating a thermionic emis- 
sion model with diffusion-limited injection currents and 
accounting for the backflow of charge at the interface. 
Such a model predicts an inverse relationship between 
the mobility and the contact resistivity, provided that the 
Schottky barrier between Au and P3HT is low. This is 
expected to be the case since the highest occupied molec- 
ular orbital (HOMO) of P3HT is estimated^ to lie be- 
tween 5.1 and 5.2 eV, close to the work function of Au 
(5.2 eV)E3- 

In the case of lower work function metals such as Cr 
and Cu (~ 4.7 eV) pj|, a significant Schottky energy bar- 
rier, A, for holes is expected to exist at the OSC-metal in- 
terface. For a channel mobility, zx, that is thermally acti- 
vated with characteristic energy E a , typically <~0.1 eV, 
the same model predicts a temperature dependence of 



the contact resistance R s cx exp([_E a + Aj/fcsT). For 
A ^ 0.3 eV, the temperature dependence of i? s is there- 
fore predicted to be much stronger than the Au elec- 
trode case. However, in recent studies of charge injection 
in both bottom-contact P3HT field effect transistors [Tc^ 
and hole injection from an Ag electrode into poly- 
dialkoxy-p-phenylene vinylene[l7j, only a weak tempera- 
ture dependence of contact resistance or injecting current 
was observed. These results imply that the diffusion- 
limited thermionic emission model is inadequate. 

In this work, we examine this issue through the tem- 
perature and field dependence of charge injection in 
bottom contact OFETs based on P3HT with different 
source/drain electrode materials. To differentiate be- 
tween the current-voltage characteristics of the channel 
and the contacts, we examine the scaling of device cur- 
rent with channel length, employ the gradual channel 
approximation^, and divide the total source-drain volt- 
age Vd into a channel component and a voltage dropped 
at the contacts, Vc ■ We assume, as supported by scan- 
ning potentiometrvp^l . that Vc is dominantly dropped 
at the injecting contact for metals with a significant A. 
We use Id — Vb data from a given series of devices of 
varying channel length, L, and fixed width, W, to ex- 
tract both [i and Id — Vc for this OSC/metal interface. 
As expected, the Id — Vq characteristics of a specific 
OSC/metal interface are unique at a given temperature 
and gate voltage, independent of L. We analyze the 
field and temperature dependence of the injected current 
through a recent analytical model 18] of charge injection 
from a metallic electrode into a random hopping system. 
With reasonable fit parameters, this model agrees well 
with the observed temperature and field dependence of 
the injected current. We also discuss the distance scale 
over which Vc is dropped, and further experimental av- 
enues to explore. 



II. EXPERIMENTAL DETAILS 

Devices are made in a bottom contact configuration^ 
on a degenerately doped p+ silicon substrate used as a 
gate. The gate dielectric is 200 nm of thermal Si02. 
Source and drain electrodes are patterned using electron 
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beam lithography in the form of an interdigitated set of 
electrodes with a systematic increase in the distance be- 
tween each pair. The channel width, W, is kept fixed for 
all devices. Three different kinds of metallic electrodes 
(Au, Cr, Cu) were then deposited by electron beam evap- 
oration followed by lift off. (25 nm of each, preceded by 
2.5 nm of Ti adhesion layer; no Ti layer for Cr samples). 
Au electrodes were cleaned for one minute in a 1:1 solu- 
tion of NH4OH: H2O2 (30%), rinsed in de-ionized water, 
and exposed for about 1 min to oxygen plasma. The Cr 
samples were cleaned in the same manner followed by a 
last step dipping in a buffered HF solution for under 10 
seconds. The HF is believed to etch the native Si02 ox- 
ide, exposing a fresh layer of dielectric. Cu electrodes on 
the other hand were only exposed to less than 25 seconds 
of O2 plasma to clean the organic residue from the lift 
off. We found that Cu samples exposed to any cleaning 
procedure except for short O2 plasma generally exhibited 
very poor transport properties. 

The organic semiconductor is 98% regio-regular 
P3HT[jij a well studied material^Ellil. As received, 
RR-P3HT is dissolved in chloroform at a 0.02% weight 
concentration, passed through PTFA 0.02 micrometer fil- 
ters, and solution cast onto the clean substrates, with 
the solvent allowed to evaporate in ambient conditions. 
The resulting films are tens of nanometers thick as deter- 
mined by atomic force microscopy. The measurements 
are performed in vacuum (~ 10~ 6 Torr) in a variable- 
temperature probe station using a semiconductor param- 
eter analyzer (HP4145B). 

III. RESULTS AND DISCUSSION 

The devices operate as standard p-type FETs in ac- 
cumulation mode[23, UM- With the source electrode 
grounded, the devices are measured in the shallow chan- 
nel regime (Vb < Vg)- Here we mainly concentrate on 
the experimental results of charge injection from Cr and 
Cu electrodes, as results on Au samples have been pub- 
lished elsewhere [5! 125). 

A. Extracting contact current-voltage characteristics 

Figures ^1 and ^> show the transport characteristics 
{I'D — Vb) of a Cr device with L = 25 and W — 200 /im 
at T = 300 K and T = 160 K for a series of V G s. The 
rise of current with voltage in the low Vb regime (Vb < 
3 V) is slightly super-quadratic. At higher drain voltages, 
the current becomes less nonlinear, at least in part due 
to saturation effects in the transistor. Our focus here 
in on low source-drain bias voltages. For comparison, 
Id — Vb plots for an Au sample with the same geometric 
parameters as above is also included in Figs.^ and^i. 

The severe nonlinearity observed in the Cr data, in 
contrast to the linear injection from Au, is traditionally 
attributed to the existence of a Schottky energy bar- 
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FIG. 1: (a) Transport characteristics (I D - Vb) of a P3HT 
OFET with Cr source and drain electrodes, with L — 25 /im, 
W = 200 fim at 300 K for several gate voltages, (b) Same 
device at 160 K. (c) Analogous plot for a sample with Au 
electrodes of the same geometry, at 300 K and (d) at 160 K. 

rier, A, between the Fermi level of the metal and the 
HOMO of the polymer. It is generally acce ptedO E3 
that the size of A alone determines the nature of trans- 
port in a device, i.e., injection limited vs. bulk (or space 
charge limited) transport. In the space-charge case, most 
relevant when A is relatively small, the easily injected 
carriers significantly alter the local electric field from 
the average field imposed by the electrode potentials, 
and correspondingly limit the current. From individual 

— Vb traces, it can be difficult to distinguish between 
space charge effects and nonlinear contact injection be- 
cause space charge limited transport is also significantly 
nonlinear |26|. One means of distinguishing the two is the 
detailed dependence of current on V and L|7lll3,l2q. A 
weak temperature dependence of current is also another 
signature of injection limited transport from metal into 
organic semiconductor!^. 

The injection properties of the contact can be 
examined^ by splitting the channel into two regimes 
of contact and the main channel. A voltage of Vc is 
dropped across the contact with the remaining V c h = 
Vb — Vc across the main channel. Using the charge con- 
trol model[2j|, Jq can be written as: 

I D = WCfx[V G -V T -V(x)}^-, (1) 

dx 

where V(x) is the potential in the channel at some posi- 
tion x, Vt is the threshold voltage, C is the capacitance 
per unit area of the gate dielectric, and fi is the intrinsic 
channel mobility. Integration of Eq. from x = to 
L — d gives: 

p^(L - d) = (V G - V T )(V D - V C ) - i(V D 2 - V c 2 ), (2) 

where Vc is dropped across d, a characteristic deple- 
tion length near the contacts. In this treatment, we 
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assume Vc to be entirely dropped across the injecting 
contact. Scanning potentiometry experiments in this 
materialjlllGIj have previously shown that, in systems 
with significant A, most of the potential drop due to con- 
tacts occurs at the source, where holes are injected into 
the channel. 

Eq. ©can be used to extract a value of Vc for any pair 
of (Vdj-Td) data, though there is no independent way of 
knowing the correct value of fi. With an array of devices, 
one can use the length dependence of Id to address this 
difficulty. At a given T and Vg, a series of Id — Vd data 
is collected from devices with different channel lengths. 
The corresponding Id — Vc is calculated from Eq. @ 
for all the different L. If the contact and channel trans- 
port properties in each device are identical, the correct 
value of fi would make all the different Id — Vc curves 
collapse onto one, since the injection characteristics of a 
particular OSC/metal interface should be unique and set 
by material properties and the (fixed) channel width and 
electrode geometry. This technique allows for the simul- 
taneous extraction of fi and Id — Vc- Since the average 
source-drain field in our devices is low (< 10 3 V/cm), 
no significant field dependence of \x is expectedp] or ob- 
served. 

To confirm this method of extracting fi and Id — Vc , 
we fabricated a series of devices (in a two-step lithogra- 
phy process) with alternating Au and Cr electrodes. The 
data is then taken twice for each device, once with the 
source electrode on Cr with the drain on Au and the sec- 
ond time vice versa. Fig. [21 shows a plot of extracted 
Id — Vc for injection from Cr and Au at T = 240 K and 
Vg = — 80 V. We noticed that there is still a minute non- 
linearity present in data for Au at that is not present in 
all-Au devices. We believe that this is consistent with 
a small contact voltage at the drain, as was seen in the 
potentiometry profile of Cr/P3HT devices in Ref. [Io| . 
The Au data in Fig. [21 have been shifted toward lower 
| Vc | by 0.5 V to account for this. The value of \i that 
collapses the different length-dependent data for injec- 
tion from Cr onto one Id — Vc curve is identical to that 
inferred from the length-dependence of the channel resis- 
tance when injection is from Au in the same devices. This 
demonstrates that this procedure of extracting Id — Vc 
is well-founded. 

The mobilities in the Au/Cr devices are lower than 
those seen in all Au or all Cr source/drain samples (dis- 
cussed below). We believe this to be due to inferior sur- 
face cleanliness of samples made in the two-step lithogra- 
phy technique. The contact resistance data for injection 
from Au agree quantitatively with the data observed in 
previous Au samples yj| . 

Fig.Elshows a plot of measured Id — Vd and the current 
corrected for contact voltages, i.e. Id — V±, for the all- 
Cr sample described above at Vc = —60 V. As seen from 
the plot, most of the total voltage is dropped across the 
contact, making these devices severely contact limited. 
For example, for a drain voltage of 2 V, Vc/Vch ~ 30. 

Fig. EJl shows the temperature dependence of \x ex- 
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FIG. 2: Extracted Id — Vc for a series of devices of width 
200 fj.m with alternating Cr and Au electrodes at 240 K and 
Vg = —80 V. Upper curve shows injection of holes from Au, 
while lower curve shows injection from Cr. Injection from Au 
is more linear and allows higher currents at lower voltages. 
The Au data have been shifted to lower |Vc| by 0.5 V to 
account for a small contact voltage at the drain. 




FIG. 3: Measured In — Vd, and the same data corrected for 
contact voltages (i.e. Id — Vch) for the Cr electrode device 
shown in Fig. [T] at 290 K and Vg = -60 V. The device is 
clearly quite contact limited. 



tracted this way as a function of T _1 for a set of devices 
with all-Cr source/drain electrodes. The temperature de- 
pendence is well approximated as thermal activation con- 
sistent with simple hopping of carriers between localized 
states in the channel. The activation energies E a for the 
mobility are quantitatively similar to those seen in all-Au 
devices The inset in Fig.^Ji shows that the activation 
energies of the injected current (at Vc = 1 V) are smaller 
than E & . In agreement with others' results [l(lll7| . this is 
inconsistent with the simple thermionic diffusion model 
of injection. As discussed in the next section, the hop- 
ping injection model predicts this weak temperature de- 
pendence. 
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FIG. 4: Temperature dependence of the extracted channel 
mobility as a function of T -1 for sets of devices with (a) 
Cr and (b) Cu source and drain electrodes at several gate 
voltages. Inset in (a): Activation energies of the injected 
current and E a as a function of Vg- Inset in (b): An example 
of "collapsed" In — Vc data for devices of different values of 
L with Cu electrodes, W = 400 fj,m, 240 K, V G = -70 V. 



A similar plot for the temperature dependence of chan- 
nel mobility in a sample with Cu source/drain electrodes 
is shown in Fig. 0Jd. The values of /i in Cu devices are 
consistently lower than those seen in Cr or Au, though 
the activation energies are very similar. The reason 
for these lower mobilities remains unclear, since con- 
tact effects have been accounted for. However, it is 
commonly observed |2t| that the field effect mobil- 
ity can be significantly different in nominally identically 
prepared samples. The values of contact resistivity, R C W 
are also higher in Cu devices. At a small Vc ~ 1 V, 
i? c W[Cu] ^3300 kil-cm, while the corresponding value 
for Cr is -RcW[Cr] ^125 kfi-cm. For devices with Au 
source/drain electrodes [5|, R C W ~ 10 kO-cm was ob- 
tained. The increase in contact resistivity from Au to 
Cu is consistent with the increase in value of A as de- 
scribed below. The inset in Fig. 2t> shows an example of 
a coalesced plot of In — Vq for different values of L for 
the Cu sample at T = 240 K and Vg = —70 V. A single 



mobility value of 3.8 x 10 3 cm 2 /Vs was used to obtain 
this collapse. 

B. Interpretation and modeling 

Having extracted nonlinear current-voltage character- 
istics for the injecting contacts, we analyze the resulting 
data in terms of a particular model of injection into disor- 
dered polymer semiconductors. As mentioned earlier, the 
diffusion thermionic models are insufficient to account 
for the weak temperature dependence of the observed in- 
jected current. The analytic treatment by Arkhipov et 
«L|l8| of charge injection from a metallic electrode into 
a random hopping system has been shown p"7ll27j to con- 
sistently explain the field and temperature dependence 
of charge injection in organic Schottky diode-like struc- 
tures. In this picture, the weak temperature dependence 
of the injection current is a consequence of a Gaussian 
distribution of statespjof. Therefore the injection pro- 
cess is easier at lower temperatures, leading to a weaker 
temperature dependence of the current. Here we apply 
Arkhipov results to our charge injection data in OFETs 
and show that this treatment, with reasonable parame- 
ters, is consistent with the measured data. 

Key to this analysis is the conclusion 0, [M| that 
a small depletion region forms in the vicinity of the 
contacts, and that Vc is dropped across this region at 
the injecting contact. We note that the values of Vc 
can be as large as a few Volts. Establishing the dis- 
tance scale, d, relevant to converting this into the elec- 
tric field at the contact is nontrivial, though reasonable 
bounds may be placed on this parameter. The resolu- 
tion of existing scanning potentiometry data in similar 
OFET structuresHEl 

establishes that d cannot exceed 
~ 400 nm. Furthermore, the lack of breakdown or ir- 
reversible device damage implies that the injecting field 
must be below the breakdown field of the OSC, so that d 
must be larger than ~10 nm. After presenting the anal- 
ysis of the Id — Vc data, we return to this issue below. 

In this Id model|l8|, the transport of carriers takes 
place in a hopping system of Gaussian energy distribution 
in close contact with the metallic electrode. This density 
of states (DOS) is given by: 

where Nt is the total spatial density of localized states, 
with a as the variance of the gaussian distribution cen- 
tered about E = 0. The emphasis is placed on the pri- 
mary injection event where a carrier from the metal is 
injected into a localized state a distance xq > a from 
the interface, where a is the intersite hopping distance. 
The potential of this carrier at any distance x from the 
interface is given by 

U(x,E) = A- — — eF x + E, (4) 

lbneoex 
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where A is the energy difference between the Fermi ievei 
of the metal and the center of DOS in the semiconductor, 
Fq is the external field at the contact, e is the elemen- 
tary charge, and e is the relative dielectric constant of 
the polymer. Once a carrier is injected into a localized 
state in the polymer, it can either go back to the metal 
due to the attractive image potential, or escape with a 
finite probability to diffuse into the bulk. The escape 
probability can be solved u sing the Id Onsager problem 
as outlined in detail in Ref. [18j . The final result predicts 
the injection current density as follows: 
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Here, v is the attempt-to-jump frequency, T is the tem- 
perature, 7 is inverse localization length, and the Boltz- 
mann function Bol(E) is defined as: 



Bol(E) = exp(-E/k B T),E > 0, 
= 1,E<0. 



To apply this model, we first need to fix the parame- 
ters (7, a, and 7. It is possible to extract a from a model 
of carrier transport in a disordered Gaussian density of 
states |32u 33] by plotting ln(/uo) vs. T -2 , where no is the 
value of the zero-field mobility. We note that our data ap- 
pear to be better described as exponential in T" 1 rather 
than T~ 2 ; nonetheless this procedure provides an esti- 
mate for a value of a. Calculation of the values of a and 
7 can be difficult, as one has to use the strong field 
dependence of mobility. As mentioned earlier, our data 
are acquired in a low enough source-drain average field 
that no field-dependence of fi may be inferred. There- 
fore, we chose a and 7 consistent with reported values in 
literature |32j or previous experiments (23. We note that 
changing a or 7 over a reasonable range mainly affects 
the overall prefactor of the current (as described below), 
without significantly altering the shapes of the predicted 
curves. 

With a, a, and 7 held fixed, the only parameters that 
can be adjusted to fit Eq. © to a plot of data are A, 
a prefactor K = AvNt (where A is an effective injec- 
tion area), and d, where Fn = Vc/d. We observed, as 
discussed in detail in Ref. [jjj , that the nonlinearity in a 
plot of Id — Vc is mainly controlled by the value of A and 
the strength of the electric field. At F ~ 5 x 10 7 V/m 
or higher, the plots are severely nonlinear and the tem- 
perature dependence of the current would be extremely 
weak. At lower fields, the nonlinearity is less severe and 
the temperature dependence is stronger. Therefore, the 
value of d is paramount, and constrained as described 
above. 




FIG. 5: Id — Vc data from a set of Cr electrode devices, 
together with curves from the model of Eq. JKJ at gate voltages 
(a) -30 V, and (b) -60 V. 



Having expected A ~ 0.3 eV, we find that d cannot be 
below 100 nm within this model; the resulting large val- 
ues of injecting field Fq would yield temperature and field 
dependences inconsistent with those observed. Table |U 
summarizes the parameters used to model the injection 
data for both all-Cr and all-Cu sets of devices. FiguresEK 
and b show plots of Jd — Vc and the corresponding nu- 
merical integration of Eq. JSJ using parameters given in 
Table [i] (with an appropriate value of prefactor K) for 
injection from Cr over a representative range of temper- 
ature for gate voltages -30 V and -60 V. Fig. |5| shows the 
temperature dependence of the injected current in low 
Vc regime and the Arkhipov fit to the data. Notice that 
the predicted temperature dependence in the diffusion 
thermionic model is much stronger than the Arkhipov 
model if the same A = 0.23 eV is used. Fig. [7\ shows a 
plot of Id — Vc for a Cu sample at Vg = —60 V. The fit 
to the data is valid only in the low Vc regime as satura- 
tion effects in the transistor start to affect Id — Vb data 
at large Vb- Note that the procedure outlined above to 
extract the Id — Vc data assumes that devices are firmly 
in the gradual channel limit, with no saturation effects. 
Also, the effects of leakage currents to the gate electrode 
in the immediate vicinity of Vc = at low temperatures 
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TABLE I: Parameters used to model the 7d — Vc data of this 
study within the charge injection treatment of Ref. (la| for 
all T. The relative dielectric constant e of the polymer was 
assumed to be 3. At each gate voltage a single numerical 
prefactor was the only necessary adjustment. 
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FIG. 6: Temperature dependence of the injected current 
for Cr electrodes at low Vc, with Eq. © fits to the data. 
The dashed line shows the predicted T dependence of the 
thermionic diffusion model for the same value of A. 

Nt is the only gate-dependent parameter in this model. 
The parameters listed in Table [I] are kept fixed for all 
temperatures and all gate voltages. Since the prefactor in 
Eq. © is the product AuNt, it is therefore not possible 
to determine an exact value for just Nt. For a cross- 
sectional area of injection of A ~ 25 nmx2 x 10 5 nm and 
v ~ 10 13 s -1 , we find N t ~ 1.1 x 10 22 cm" 3 for V G = 
-60 V and N t ~ 8.3 x 10 21 cm" 3 for V G = -30 V. These 
values are consistent with other experiments [Til l2?| . 

Tablc[I]shows that the obtained injection barrier height 
for copper is about 80 meV higher than that for Cr. This 
difference in the barrier energy is not unreasonable, and 
may be attributed|lO| to an interfacial dipole layer at the 
interface changing A by a small amount. In-situ ultra- 
violet photoemission spectroscopy measurements would 
be well-suited to testing this hypothesis |34l l35l l36 | . The 
higher injection barrier for copper is consistent with the 
observed higher contact resistivity and lower overall cur- 
rents observed in Cu. 

The length of the presumed depletion region is also a 
bit higher in Cu samples (by ~ 80 nm), though it does 
not necessarily reveal why the mobility is lower in these 
devices. The origin of these depletion regions in the vicin- 
ity of the contacts is not understood in detail. Recent 
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FIG. 7: Id — Vc for devices with Cu electrodes (W = 400 nm) 
at Vg = —60 V over a representative temperature range, to- 
gether with results from Eq. and the parameters of Tabled] 



2d electrostatic modeling [3l| of OFETs has shown that 
the effect of significant energy barriers at the injecting 
electrode is formation of regions of low carrier concen- 
tration (and mobility) near the contacts. These studies 
place the extent of these regions at about 100 nm from 
the contact, depending on Vg- Another possible origin 
for regions of reduced mobility near metal contacts with 
significant barriers is charge transfer and band bending 
near the interface. Since conduction in these materials 
can be treated as percolative variable range hopping 37], 
(i is a natural function of the density of available hop- 
ping sites. The occupation of those sites can be strongly 
modified by interfacial charge transfer between the metal 
and the OSC. Improved local probes (nm-resolution scan- 
ning potentiometry, cross-sectional scanning tunneling 
microscopy) would be extremely useful in better under- 
standing these depletion regions. 



V. CONCLUSIONS 

Transport properties of a series of organic field effect 
transistors with P3HT as the active polymer layer and 
Cr, Cu and Au as the source/drain electrodes were ex- 
amined over a temperature range. The contact current- 
voltage characteristics for these devices were extracted 
from the length dependence of conductance, with the as- 
sumption that the injection barrier primarily applies to 
holes being injected from the source. This procedure was 
checked for consistency using devices with electrodes of 
alternating metal composition. The data confirm that 
the weak temperature dependence of the injected current 
cannot be simply explained using the general diffusion- 
thermionic emission models. With reasonable values of 
parameters, a more sophisticated model of hopping in- 
jection into a disordered density of localized states, with 
emphasis on the primary injection event, is consistent 
with the field and the temperature dependence of the 
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data over a broad range of temperatures and gate volt- tion. 
ages. 
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